Numerical analysis of the planewave discretization of orbital-free and Kohn-Sham models Part I: The Thomas-Fermi-von Weizsäcker model
نویسندگان
چکیده
We provide a priori error estimates for the spectral and pseudospectral Fourier (also called planewave) discretizations of the periodic Thomas-Fermivon Weizsäcker (TFW) model and of the Kohn-Sham model, within the local density approximation (LDA). These models allow to compute approximations of the ground state energy and density of molecular systems in the condensed phase. The TFW model is stricly convex with respect to the electronic density, and allows for a comprehensive analysis (Part I). This is not the case for the Kohn-Sham LDA model, for which the uniqueness of the ground state electronic density is not guaranteed. Under a coercivity assumption on the second order optimality condition, we prove in Part II that for large enough energy cut-offs, the discretized Kohn-Sham LDA problem has a minimizer in the vicinity of any Kohn-Sham ground state, and that this minimizer is unique up to unitary transform. We then derive optimal a priori error estimates for both the spectral and the pseudospectral discretization methods.
منابع مشابه
Numerical analysis of the planewave discretization of some orbital-free and Kohn-Sham models
We provide a priori error estimates for the spectral and pseudospectral Fourier (also called planewave) discretizations of the periodic ThomasFermi-von Weizsäcker (TFW) model and for the spectral discretization of the Kohn-Sham model, within the local density approximation (LDA). These models allow to compute approximations of the ground state energy and density of molecular systems in the cond...
متن کاملMagnetic Thomas-Fermi-Weizsäcker model for quantum dots: a comparison with Kohn-Sham ground states
The magnetic extension of the Thomas-Fermi-Weizsäcker kinetic energy is used within densityfunctional-theory to numerically obtain the ground state densities and energies of two-dimensional quantum dots. The results are thoroughly compared with the microscopic Kohn-Sham ones in order to assess the validity of the semiclassical method. Circular as well as deformed systems are considered. PACS. 7...
متن کاملA perturbation-method-based post-processing for the planewave discretization of Kohn-Sham models
In this article, we propose a post-processing of the planewave solution of the Kohn–Sham LDA model with pseudopotentials. This post-processing is based upon the fact that the exact solution can be interpreted as a perturbation of the approximate solution, allowing us to compute corrections for both the eigenfunctions and the eigenvalues of the problem in order to increase the accuracy. Indeed, ...
متن کاملConjugate-gradient optimization method for orbital-free density functional calculations.
Orbital-free density functional theory as an extension of traditional Thomas-Fermi theory has attracted a lot of interest in the past decade because of developments in both more accurate kinetic energy functionals and highly efficient numerical methodology. In this paper, we developed a conjugate-gradient method for the numerical solution of spin-dependent extended Thomas-Fermi equation by inco...
متن کاملFinite-temperature orbital-free DFT molecular dynamics: Coupling Profess and Quantum Espresso
Implementation of orbital-free free-energy functionals in the Profess code and the coupling of Profess with the Quantum Espresso code are described. The combination enables orbital-free DFT to drive ab initio molecular dynamics simulations on the same footing (algorithms, thermostats, convergence parameters, etc.) as for Kohn–Sham (KS) DFT. All the non-interacting free-energy functionals implem...
متن کامل